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Two Functions Whose Powers Make Fractals 


Marc Frantz 


1. AN OLD FRIEND /. In [5] Richard Darst and Gerald Taylor investigated the 
differentiability of functions f p (which for our purposes we will restrict to (0, 1)) 
defined for each p > 1 by 

(0 if x is irrational 
J \ x ) | l/n p if x = m/n with (ra, n) = 1. 

These functions are powers of the function f = f 1 , which frequently appears in 
real analysis courses because of its interesting property of being discontinuous at 
rational points and continuous at irrational points. Darst and Taylor showed that if 
1 < p < 2, then f p is nowhere differentiable, and if p > 2, then f p is differen- 
tiable almost everywhere. Graphs of / 1 and its variants (Figure 1) suggest there 
may be fractals involved somewhere, and indeed this is the case. 




(b) 


Figure 1 . (a) The graph of f 1 . (b) A graph of f p for p > 1. 


In this article, “fractal” means a subset of U n for some n , whose Hausdorff 
dimension is not an integer. For a discussion of Hausdorff measure and dimension 
see [6, Chapter 1]. Although the graphs of the functions f p are not themselves 
fractals in this sense, they do display an interesting kind of self-similarity. To be 
precise, if G p is the graph of f p for some p > 1, and the functions L p , R p , and T p 
are defined on [0, 1] X [0, 1] by 


L p (x,y) 


IJL 2L_) 

\ i +*’ (i +x) p )’ 


R P ( x >y) 


1 +x’ (1 + x) p I’ 


T p (x,y) = ( 1/2,1/20, 

then L p (G p ), R p (G p ), and T p (G p ) are the left half, right half, and top point, 
respectively, of G p . This is left as an exercise to the reader. In fact, the functions 
L p , R p , and T p can be used to sketch the graphs G p by an iterative procedure 
essentially identical to the iterated function system method made popular by 
Michael Barnsley in [1]. For more information on graphing, see the remarks at the 
end of this article. 
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The real fractals show up when we deal with differentiability. Fortunately, we 
can avoid the usual difficult estimates for Hausdorff dimension by employing 
theorems from the theory of rational approximation of real numbers; this theory is 
closely related to the differentiability of the functions f p . In the language of that 
theory, an irrational number x is said to be approximable by rationals to order p 
[7, p. 88] if there exists a positive constant c, depending only on x, such that the 
inequality 

\x - m/n\ < cn' p (1) 

has infinitely many rational solutions m/n with n > 0 and ( m,n ) = 1. Assuming 
that x,m/n e (0, 1), we have f p (x) = 0 and f p (m/n ) = n~ p , so condition (1) is 
equivalent to the existence of a positive constant c, depending only on x, such that 


f p O) - f p (m/n ) 


x — m/n 


n~ p 1 

\x — m/n\ c 


( 2 ) 


for infinitely many rationals m/n , arbitrarily close to x. But x, being irrational, is 
a minimum point of f p , so (f p )'(x) exists if and only if (f p )'(x) = 0. Hence (2) 
implies that (f p )'(x) does not exist, and the reverse implication is clearly true as 
well. We have proved 


Lemma 1. An irrational number x <e (0, 1) is approximable by rationals to order p if 
and only iff p is not differentiable at x. 

We refer to the set of points at which f p is nondifferentiable as Nondiff f p , and 
we denote the Hausdorff dimension of a set A by dim A. 


Theorem 1. For p > 2, dim Nondiff f p = 2/p. 


Proof: It suffices to prove the theorem for the set N p of irrationals in Nondiff f p . 
To show that dim N p > 2/p , we need only exhibit a subset A of N p such that 
dim A = 2/p. The existence of such a set was established by Jarnik ([6, Theorem 
8.16(a)], [8]), who proved that if p > 2 and A is the set of irrational numbers x 
such that 

\x — m/n\ < n~ p 

for infinitely many rationals m/n , then dim A = 2/p. By Lemma 1, f p is not 
differentiable at any such x, so A c N p . 

To show that dim N p < 2/p , fix an arbitrary s with 0 < e < p — 2. If x e N p , 
then there must exist a positive constant c = c(x ) such that for infinitely many 
rationals m/n , 

\x — m/n \ < c/n p < n e /n p , 


the latter inequality being true for sufficiently large n. Thus for all x ^ N p , we 
have \x — m/n\ < n~ {p ' e) for infinitely many n, so Jarnik’s theorem implies that 
dim N p < 2/(p - s ). Since s was arbitrary, it follows that dim N p < 2/p. m 

It is hopeless to look for second derivatives of the functions f p , but we can 
come close. For sufficiently large p , the related limit 


/£(» = lim 

h -»0 


fP(x + h) - 2fP(x) +fP(x-h ) 

h 2 


( 3 ) 


exists, and fractals are again involved. 
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Theorem 2. Let N* c (0, 1) be the set of points x for which f%(x) does not exist. 
Then for p > 2, N* = Nondiff f p/1 . Thus , by Theorem 1, dim N* = 4/p for p > 4. 

Proof: It should be clear that f%(x) exists if and only if x is irrational and 
/*(•*:) = 0. Furthermore, since f p (x) = 0 for every irrational x, the quotient in (3) 
is nonzero only if one of ^ + h or jc - h is rational, and it must be exactly one of 
these, for if both were rational, then x would be rational, being the average of two 
rational numbers. Hence the quotient in (3) is nonzero if and only if exactly one of 
x + h or x — h \s equal to m/n for (m, n) = 1, and then 
f p (x + h)-2f p (x) +f p (x-h) = l/n p 

h 2 ( x - m/n) 2 

= l 1 /n p/1 \ 2 ( f p/2 (x) ~f p/2 (m/n) \ 2 

\x — m/n J \ x — m/n J 
Thus f%(x) fails to exist precisely when x e Nondiff f p/ 2 . ■ 

2. A RELATED APPLICATION. When pondering counterexamples such as f 1 
and its powers, it is natural to wonder if such things might apply to the real world. 
In fact, graphs like those in Figure 1 appeared some time ago in an applied setting 
(Figure 2). 



Figure 2. The function AO. After Jensen, Bak, and Bohr [10]. 


Jensen, Bak, and Bohr [10] made use of the graph shown in Figure 2 (cf. [10, 
Fig. 5]) to help illustrate the behavior of an important dynamical system — the 
so-called critical circle map — which is used in the analysis of many physical 
systems, ranging from the driven, damped pendulum to phenomena in condensed- 
matter physics (cf. [10]). Figure 2 shows the left half of the graph of a function 
Afl : [0, 1] -> U that was introduced in [10]. Notice that this graph bears a striking 
resemblance to the left half of Figure 1(b). In both cases the right half is the mirror 
image of the left. The interpretation of the graph is roughly as follows. The 
independent variable (the winding number) describes the periodic or nonperi- 
odic behavior of the system, and x in turn is determined by a parameter f 1 (not 
depicted in the figure), which can represent some physical property of the system. 
If this parameter changes, what happens to xl That depends on AIK*). If, say, 
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A ft(x) = 0.12, then ft can vary over an interval of length 0.12 without affecting x; 
the system is then called “ mode-locked”. On the other hand, if Af l(x) = 0, then 
the slightest variation in ft will cause x to change to a new value, and the behavior 
of the system changes. No exact formula is known for Aft(x), but it is known that 
Aft(x) = 0 for each irrational x e [0,1] and Af 1(m/n) > 0 for each rational 
m/n [0,1]. Moreover, numerical evidence presented in [10] suggests that 0 < 
Aft(x) < Af Pl (x) for all x e (0, 1), where A is a positive constant and p 1 > 2. If 
this is true, then the result of Darst and Taylor — as well as Theorem 1 — imply that 
Aft is differentiable almost everywhere on [0, 1]. Although it is known that Aft 
cannot be equal* to Af p for any constants A and p, further numerical results in 
[10] suggest that Af 229 , for some constant A, could serve as a rough model for 
Aft. In view of the tractability of the functions f p , and in view of interest in the 
Hausdorff dimension of sets related to Aft, this could be a promising avenue of 
investigation; see [2], [3], [9], [10]. 

The function Aft is, in fact, a step-width function for a devil’s staircase 
associated with the critical circle map (Figure 3). The function W : [0, 1] -> R, 


0.8 h 
0.6 / 

W{Vt) r . 

0.4 

0.2 ■ “ 

o.o ■ 

0.0 0.2 0.4 0.6 0.8 

n 

Figure 3. Mode-locking staircase. After Jensen, Bak, and 

whose staircase graph appears in the figure, takes the previously-mentioned 
parameter ft as its input; the output IT(ft) is the winding number associated with 
ft. The value Aft(x) gives the width of the step on the staircase with ordinate 
x = IT(ft); that is, Aft(x) = A(IT _1 ({x})), where A denotes Lebesgue measure. 

Results in [ 10 ] indicate that the Hausdorff dimension of the corresponding 
Cantor set (the closure of the set of ft-values with irrational winding numbers) is 
approximately 2/2.29, which by Theorem 1 is the same as dimNondiff f 229 . Is this 
a coincidence? Probably not. In the next section we use the step-width function 
approach with respect to a class of Cantor functions studied recently by Darst, and 
show that each corresponding Cantor set has the same Hausdorff dimension as the 
nondifferentiability set of the associated step-width function. 

3. A RELATED FUNCTION g WHOSE POWERS CORRESPOND TO CANTOR 
SETS. In [4] Darst describes a class of Cantor sets C a as follows. For each 
a e (0, 1/2), put b = 1 - 2a and delete an open interval of length b from the 
center of [0, 1], leaving two closed intervals of length a. From the center of each of 


1.0 


Bohr [10]. 


*For fixed n , the values Af p (m/n ) are the same whenever ( m,n ) = 1. The corresponding values 
A Cl(m/n) are not all equal; it is their average that is well-approximated by Af 229 (m/n). See [10] for 
details. 
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1/4 1/2 3/4 1 0 i/4 1/2 3/4 1 

. . . drop the “steps” ... ... mark the tops. 


Figure 4. Obtaining the step-width function h 1/2t from the classical Cantor function <p 1/3 . For any 
a e (0, 1/2), the corresponding Cantor set C a satisfies dimC^ = dimNondiff h a . 


these intervals, delete a scaled version of the first open interval; that is, delete an 
interval of length ab , leaving four closed intervals of length a 1 . Repeating the 
process indefinitely leaves us with the set C a \ in particular, C 1/3 is the classical 
Cantor middle-thirds set. Darst then defines corresponding Cantor functions 
cp a : [0, 1] -> [0, 1] in the obvious way. The graphs of these functions, including that 
of the classical Cantor function <p 1/3 , are staircases with steps of finite width at 
each ordinate m /2 n e (0, 1) where m is an odd integer, and steps of zero width 
elsewhere. For our purposes, let us ignore the “steps” at ordinates 0 and 1, and for 
each a e (0,1/2) define a step-width function h a : (0, 1) -> U by h a (x) = 
hicp^dx})) (see Figure 4). Where the steps have finite widths, they are the lengths 
a n b of deleted intervals in the construction of C a . A little thought shows that 
h a (x ) = a n ~ 1 b if x = m/2 n with m odd, and h a (x ) = 0 otherwise. If we set 
p a = ln(l/a)/ln2 = log 2 (1/a), we can write 


h 


a 


(b/a)( l/2 n ) Pa if x = m/2 n with m odd 
0 otherwise . 


Now the rationals m /2 n are sometimes referred to as the dyadic rationals. If, in 
our definition of the functions f p , we replace the rationals with the dyadic 
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rationals, we get analogous functions g p : (0, 1) -> R defined by 

(l/2 n ) p if x = m/2” with m odd 
0 otherwise . 

for each p > 1. Analogously, the function g = g 1 and its powers are discontinuous 
at and only at each dyadic rational point of (0, 1). Notice that h a = ( b/a)g Pa , and 
as a runs over (0, 1 /2), the values p a run over (1, oo). Thus the functions g p for 
p > 1 correspond in a precise one-to-one way to the Cantor sets C a . 

In [4] Darst addresses the nondifferentiability of the Cantor functions cp a . He 
defines their nondifferentiability points as those points where there is no deriva- 
tive, finite or infinite, and shows that the corresponding set has Hausdorff 
dimension [ln(2)/ln(l/^)] 2 . He also shows that for each a e (0,1/2), dim C a = 
ln(2)/ln(l/fl). The next result fits nicely with these. 

Theorem 3. For each a e (0,1/2), dim Nondiff = ln(2)/ln(l/«). Equivalently , 
for each p > 1, dimNondiff g p = 1/p. 

Proof: Clearly Nondiff h a = Nondiff g p % and l/p a = ln(2)/ln(l/a). Thus it suf- 
fices to show that dim Nondiff g p = 1/p for each p > 1. Moreover, following the 
proof of Theorem 1, it suffices to prove this result for the subset M p of points in 
Nondiff g p that are not dyadic rationals. To show that dim M p > 1 /p, we need 
only exhibit a subset B of M p such that dim B = 1 /p. The existence of such a set 
follows from a result by Falconer [6, Theorem 8.15], which implies that if p > 1 
and B is the set of numbers x that are not dyadic rationals and satisfy 

\x — m /2" | < (2”)~ p (4) 

for infinitely many dyadic rationals m/2”, then dim B = 1/p. To see how this 
applies, observe that x is a minimum point of g p , so (g^OXx) can exist only if it is 
zero. But inequality (4) implies 

gPjx) -g>(m/2 n ) (2' I )~ P ; 1 

x - m/2” |jc - m/2” | - 

for infinitely many dyadic rationals m/2”. Thus g p is not differentiable at any 
ig 5, and B c M p . 

To show that dim M p < 1/p, fix an arbitrary e with 0 < s < p — 1. If r e M p , 
then there must exist a positive constant c = c(x) such that for infinitely many 
dyadic rationals m/2”, 

\x - m/2” | < c/(2") p < (2”)V(2”) p , 

the latter inequality being true for sufficiently large n. Thus for all x ^ M p , we 
have \x - m/2”| < (2”) _(/?_£) for infinitely many n, so Falconer’s result implies 
that dim M p < l/(p — s). Since s was arbitrary, it follows that dim M p < 1/p. 

■ 

The functions g p make fractals, except for g 1 : Nondiff g 1 is all of (0,1). 
However, the graph of g 1 is very fractal-like (Figure 5). If you think you see a right 
isosceles Sierpinski triangle there, it’s because the points of the graph above the 
x-axis are precisely the apexes of those subtriangles of the Sierpinski triangle 
whose hypotenuses lie on the x-axis. The picture can also be used to give an easy 
proof of the nondifferentiability of g 1 ; one can actually see its nondifferentiability. 
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Theorem 4. The function g 1 is not differentiable at any point of (0, 1). 


Proof: Since g 1 is not continuous at any dyadic rational, let us focus on the other 
points. Choose any point x e (0, 1) that is not a dyadic rational. As in the previous 
proof, (g 1 )/.*) exists only if (g 1 )/*) = 0; we show that this cannot happen. Draw 
the triangle Aafiy, where a = (x, 0), /3 = (0, 1/2), and y = (1, 1/2). Thinking of 
the top as the base, A a /3y has base 1 and height 1/2, so at any ordinate 1/2", the 
triangle has width 2/2". Additionally, at any ordinate 1/2", the graph of g 1 has 
points with abscissas 1/2", 3/2", . . . , (2" - l)/2", which are spaced 2/2" units 
apart. Thus at every ordinate 1/2", triangle Aafiy contains a point of the graph. 
Since its slanted sides have slope at least 1/2, it follows that (g 1 )'(jc) cannot be 
zero. ■ 


Remarks. An avenue for future work is to think of the functions f p for /? > 2 as 
step-width functions, construct corresponding Cantor sets K p and Cantor func- 
tions ifj p , and show (presumably) that dim K p = 2/p. Staircase graphs of the 
functions p p would of course have steps of finite width at every rational ordinate, 
making them similar to the mode-locking staircase in Figure 3. Put differently, we 
can start with f p for p > 2, and in effect reverse the process in Figure 4. This is 
possible because f p (m/n ) = l/n p for precisely (pin ) rationals m/n e (0, 1), where 
<p denotes the Euler totient function. Hence the total length c p = E co n = l pin)/n p 
of the steps is finite, since p > 2 and pin) < n for each n. If desired, the total 
length can be scaled to 1 by using the functions (1 /c p )f p as step-width functions, 
just as the functions h a were defined as suitable multiples of the functions g Pa . 

Since the sets Nondiff f p and Nondiff g p , for sufficiently large p, are fractals by 
the definition used here, it is natural to wonder whether they might be indicated 
graphically in some way, as is often done with the Cantor middle-thirds set. 
Unfortunately, even without the rationals, both sets are dense in (0, 1), so no 
helpful illustrations are possible. To prove this for Nondiff g p , one can start with 
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an arbitrary open interval in (0, 1), and construct in it a nondifferentiable irrational 
point x by requiring x to have a binary “decimal” expansion which ultimately 
consists of longer and longer strings of 0’s separated by single l’s. The growth of 
the strings must depend appropriately on p; the details are left to the interested 
reader. 

The density of the irrational subsets N of the sets Nondiff f p has an interest- 
ing reason behind it. Recall that the Liouville numbers are irrational numbers x 
such that for each positive integer k , there exist integers m and n with n > 1 such 
that 


\x — m/n\ < n k . 

Since we may assume that ( m , n) = 1, it is easy to see that x is not a differentiable 
point of f p for any p. Consequently, the Liouville numbers form a subset of 
Op > 2 Np. Conversely, if x e f] p > 2 then for each positive integer k we have 
x e N 2+k , so there must exist a positive constant c = c(x, k) such that 

\x - m/n\ < cn~( 2+k) = ( cn~ 2 )n~ k 

for infinitely many rationals m/n. Thus \x — m/n\ <n~ k for some rational m/n , 
and hence x is a Liouville number. It follows that 0 P > 2 N is the set of Liouville 
numbers in (0, 1). The Liouville numbers are transcendental numbers that form a 
first category subset of R, and there are uncountably many of them in every open 
interval; thus each set N p is also dense in (0, 1). These facts are proved in [11, 
Chapter 2], where it is also shown that the set of Liouville numbers has Hausdorff 
dimension zero. This last fact is easily deduced from Theorem 1, since 
dim Nondiff f p [ 0 as p T oo. 

Finally, a word is in order concerning the graphs of the functions f p and g p . 
Many texts include / 1 as an example, but the vast majority do not provide a 
detailed graph (an exception being [12, p. 129]). Perhaps it is thought that such 
graphs would exhibit no interesting structure, or that they would be hard to 
generate. Neither of these assumptions is true. The graphs of the functions f p and 
g p can be used to link important basic concepts of real analysis with notions of 
rational approximation, self-similarity, fractal geometry, and dynamical systems, 
and hence can serve to motivate students in the classroom. As for generating the 
graphs, there are several easy ways to do it. I would be happy to send the 
appropriate Maple or Mathematica code in response to email requests. 

ACKNOWLEDGMENTS. I am indebted to David Stern, who first pointed out to me the nondifferen- 
tiability of f l , and to Yuri Abramovich for his encouragement and advice, including the suggestion that 
fractals might lurk somewhere. 
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